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Abstract. We consider a family of non-convex integral functionals 

\ [ W{x/e,ld+hVg{x)) Ax, g G W 1 ' p (n-,R n ) 
h 2 Jn 

where W is a Caratheodory function periodic in its first, and non-degenerate 
in its second variable. We prove under suitable conditions that the r-limits 
corresponding to linearization (h — > 0) and homogenization (e — > 0) commute, 
provided W is minimal at the identity and admits a quadratic Taylor expansion 
at the identity. Moreover, we show that the homogenized integrand, which is 
determined by a multi-cell homogenization formula, has a quadratic Taylor 
expansion with a quadratic term that is given by the homogenization of the 
second variation of W. 
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1. Introduction 
In this contribution we consider periodic integral functionals of the type 

(1.1) J W(x/e,Vu(x)) dx, ttE^njR") 

n 

where e is a positive parameter, 51 an open, bounded Lipschitz domain in M 71 with 
n > 2 and W : E n xM" ->■ [0, +oo) is a Caratheodory function Y := [0, l)"-pcriodic 
in its first variable. 

Functionals of this type model manifold situations in physics and engineering. 
We are particularly interested in applications to elasticity. In this context the in- 
tegral in (jl.ip is the elastic energy of a periodic composite material with period e 
deformed by the map u : VI — > R™. In situations where the deformation is close to a 
rigid deformation, say |Vu — Id| ~ h, it is natural to introduce the scaled displace- 
ment g(x) := h^ 1 (u(x) — x) and to consider the scaled, but equivalent energy 

l e ' h {g) -=^2 J W(x/e,ld+hVg(x)) dx, g € W^ p ((l;R n ). 
n 

For small parameters e and h we expect that the energy T £,h can be replaced by 
an effective model that is simpler than the original one, but nevertheless captures 
the essential behavior from a macroscopic perspective. In this context the limit 
h — > corresponds to linearization, while the theory of homogenization studies the 



Date: May 17, 2011. 



2 



S. MULLER & S. NEUKAMM 



asymptotic behavior of the energy as e — > 0. Obviously, there are two orderings to 
derive an effective model by consecutively passing to the limits e — > and h — > 0. 
The goal of this paper is to show that both ways lead to the same limiting energy. 
In other words we prove that homogenization and linearization commute. 

It is well-known that in the case where the integrand W is convex in its second 
variable the homogenization of (jl.l|) is the integral functional Jq W-^^CVu^x)) dx, 
where the homogenized integrand is given by the one-cell homogenization formula 

WHI(F) := inf | J W(y, F + V^(y)) dy : <p £ W^(Y 

This result goes back to P. Marcellini [15, and was extensively studied with various 
methods (cf. e.g. [8], [2]). The more general setting of monotone operators and non- 
periodic dependency on the spatial variable was studied by L. Tartar [20l ITTT [2T] , 
On the other side, for non-convex W it turns out that in general the relaxation of 
W over one — or even an ensemble of finitely many copies of the periodicity cell 
Y is not sufficient for homogenization. Nevertheless, A. Braides [H [5] and the first 
author [18| showed in the 1980s that in the non-convex case the homogenization of 
the functional in leads to the multi-cell homogenization formula 




W£ZHF) := mf inf { £ / W(y 7 F+ V<p(y)) dy : <p € W^(kY;R n ) 



K kY 



We state this result in a precise manner in Theorem 16.11 below. The fact that the 
one-cell formula is not sufficient in nonlinear elasticity has been observed earlier in 
the engineering literature (see [UH2]). These authors observed that instabilities 
arise even when the one-cell formula predicts no loss of stability. Later in |13j this 
has been systematically investigated. 

In this contribution we are interested in integrands of the following type: For 
a > and p G (l,oo) let W(a,p) denote the class of Caratheodory functions 
W : R"xM" — > [0, oo) that are F-periodic in the first variable and that satisfy 
the following conditions (Wl) - (W4). We suppose that W is locally Lipschitz 
continuous and satisfies standard growth conditions of order p: 

j ±\F\ P -a<W(y,F) < a(l + \F\ P ) and 

I \W(y, F) - W(y, G)\ < o(l + \Ff- 1 + {Of- 1 ) \F - G\ 



(Wl) 



for all F,G £ M. n and a.e. y € R™. We suppose that F = Id is a natural state and 
W is non-degenerate, i.e. 

(W2) W(y, Id) = for a.e. y e R n , 

(W3) W(y, F)>\ dist 2 (F, SO{n)) for all F € M" and a.e. y £ R". 

Furthermore, we suppose that W admits a quadratic Taylor expansion at Id in the 
sense that 

/ttti \ n nc n v \W(y,Id+G)-Q(y,G)\ 
(W4) 3Q € Q : hm sup ess sup 3 = 

G^O ugf \G\ 
G^O 1 1 

Here and below Q denotes the set of Caratheodory functions 
Q : R"xM" — > [0, +oo) that are F-periodic in the first variable, non-negative 
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and quadratic in the second variable and bounded in the sense that 
VG e M" : ess sup Q(y,G) < c\G\ 2 

for a suitable constant c > 0. 

Remarks. (1) In elasticity it is natural to assume that W is frame indifferent, 
i.e. 

W(y, RF) = W(y, F) for all R e SO(n), F £ M" and a.e. y e R n . 

We have not introduced this assumption explicitly, because it is not required 
in the proof. Nevertheless, the lower bound (W3) is motivated by frame 
indifference. 

(2) (W2) & (W3) are the main assumptions in our analysis. For energy den- 
sities representing composite materials, condition (W2) requires that each 
of the single components of the composite is stress free in a common refer- 
ence configuration. Note that this rules out the application to prcstressed 
composites, as considered in Section [7721 

(3) The combination of (W3) & (W4) can be interpreted as a generalization of 
Hooke's law to the geometrically nonlinear setting, in the sense that both 
conditions together imply that for infinitesimal small strains the stress is 
proportional to the strain. 

For W G W(a,p) we study the behavior of the homogenized integrand 
near the identity. Our first main result is the following: 

Theorem 1.1. Let W € W(a,p). Then 



(1.2) limsup 



wiZ\^+G)-Qi 1 l(G) 



\G\ 2 



0. 



The statement remains valid if we drop condition (Wl). 

This means that whenever the integrand W admits a quadratic expansion at 
Id with a quadratic term Q, then also the homogenized integrand W^™^ has a 
quadratic expansion where the quadratic term is given by the homogenization of Q. 
Partial results in this direction under strong implicit assumptions on the minimizers 
for the cell problems have been obtained in [T3"] . 

Remark 1.2. Conditions (W3) & (W4) imply that the quadratic integrand Q{y,-) 
is a positive semi-definite quadratic form for almost every y G K™. More precisely, 
the non-degeneracy condition (W3) implies that Q(y, •) restricted to the subspace 
of symmetric nxn matrices is positive definite; in particular, we have 

Q(y, G) > i |symG| 2 for all G € M™ and a.e. y e K". 

In virtue of Korn's inequality this guarantees that the minimization problem 

J Q(y,G + V(p(y)) dy subject to Jtpdy = 

Y Y 

has a unique minimizer in W^ 2 (Y; R n ). If the energy density W additionally satis- 
fies (W2) and is frame indifferent, then the associated quadratic form Q vanishes for 
skew symmetric matrices. In this case Q and its homogenization Qhom are energy 
densities as they typically appear in linear elasticity. 
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Theorem 11.11 shows that the expansion property (W4) is stable under homog- 
enization for energy densities in W(a,p). The next result states that also the 
non-degeneracy condition is stable under homogenization: 

Lemma 1.3. Let W <E W(a,p), then W^om e VV(o',p) for a positive constant 
a' = a'(a,n). Additionally, the map Wf!^ : M" — > [0,oo) is continuous and 
quasiconvex. 

In the language of T-convergence Theorem O implies that the T-limits of (l e ' h ) 
corresponding to linearization and homogenization commute. In order to state this 
in a precise manner, let 7 be a closed subset of dfl with positive n— 1-dimcnsional 
Hausdorff-measure. We denote the space of scaled displacements that satisfy the 
Dirichlet boundary condition on 7 by 



A := {g£ W 1,a (fi;K B ) : 5 = 0117}. 



For simplicity we assume that 7 is regular enough to guarantee that A^DW '°° (CI; R") 
is strongly dense in Ay. We consider the following functionals from L 2 (f2;R") to 
[0,+oc]: 



2°(ff) 



J W(x/e, Id +h Vg(x)) dx if g € A 

+ 00 else, 

£2 / ^om (W +h Vg(x))dx if g e A, 

n 

+ 00 else, 
y Q(x/e, Vg(x)) da; if 5 € Ay 



12 

+ 00 



else, 

Q£L(Vfl(s)) dz ifffeAy 



12 

+ CO 



else. 



Our second main result is the following: 

Theorem 1.4. LetW € W(a,p) withp > 2. TTien the following diagram commutes 



(1.3) 





(1) 


-> ?L 


(2) l 






(3) 


jh 
-^hom 


(4) 


-> 2 






where (1), (4) and (2), (3) mean Y -convergence with respect to strong convergence 
in L 2 (n:M. n ) as h — > and e — >• 0, respectively. 

In the diagram the T- limit (1), which corresponds to linearization of a hetero- 
geneous energy, was treated by G. Dal Maso, M. Negri and D. Percivale in [9]. In 
Section[5]we give a slight variant of their argument, which is adapted to assumption 
(W4). 
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Remark 1.5. On the level of the energy functional X E > , it is natural to study 
also the T-convergence behavior as e and h simultaneously converge to 0. This 
corresponds to a "diagonal limit" in diagram (11.31) . Indeed, in [T7] the second 
author proves (based on two-scale convergence methods) that I e,h F-converges to 
1° as (e, h) — > (0, 0) simultaneously. Theorem ll.4l is also related to recent works by 
A. Braides and B. Schmidt [6, 19] where the passage from pair-interaction atomistic 
models to linear elasticity is studied. In contrast to the setting considered in the 
present paper, where e describes the length scale of the material heterogeneity, 
in [5J HH] the small scale e originates from the discrete nature of the atomistic 
model and measures the interatomic distance. In |19j the passage from-discrete- 
to-continuous is obtained (in the regime e -C h 2 ) as a "diagonal limit", i.e. as a 
T-convergence result as e and h simultaneously converge to 0. It is observed that 
the derived r-limit coincides with the model obtained by linearizing the continuum 
model derived from the atomistic energy by applying the Cauchy-Born rule. The 
central assumption is a discrete version of the non-degeneracy condition (W3). 
In the atomistic setting this ensures that for sufficiently small displacements no 
oscillations on the length scale of the lattice emerge. 

The F-convergence result is complemented by the following equi-coercivity state- 
ment: 

Proposition 1.6. Suppose that W £ W(a,p) with p > 2. Then there exists a 
positive constant c\ such that 

min { l s ' h (g), lf in (g), l£ oni (g), l°(g) } > c^(g). 

for all e, h > and g 6 L 2 (Q; R"), where 

IMIw>.'(fi;R«) tfg^A-f 
+ oo else. 



Ha) 



Remark 1.7. Obviously, the map ^ (and the restriction ) is coercive and lower 

semicontinuous in L 2 (fl;M. n ) (and coercive and lower semicontinuous with respect 
to weak convergence in W 1,2 (Q,\ W 1 ) respectively); thus, the previous proposition 
implies that the energies X s ' , X^ n , X^ om and 1° are equi-coercive in the strong 
topology of L 2 {n;R n ). 

Remark 1.8. The local p-Lipschitz condition in (Wl) can be dropped (see M. Bai'a 
and I. Fonseca [3]). Nevertheless, it is a natural property of quasiconvex or even 
rank-one convex functions: Let p g [l,oo) and / : M" — > [0, oo) with f{F) < 
C(l + \F \ p ). If / is rank-one convex, then we have 

\f(F) - f(G)\ < C(l + {Ff- 1 + IGf- 1 ) \F-G\. 

We would like to remark that the p-growth and coercivity condition is too restric- 
tive for a direct application to finite elasticity, because we expect the behavior 
W(y, F) —> +oo as F for realistic materials. In this direction, in [17] the sec- 
ond author considers periodic integrands W that only satisfy (W2) - (W4). It turns 
out that Theorem 11.41 remains valid, when X^ om is replaced by the lower r-limit of 
(l e ' h ) e as e -> 0. 

The plan of the paper is as follows. In Section [2] we prove Lemma 11.31 and 
show that the non-degeneracy condition (W3) is stable under homogenization. In 
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Section [3] we prove Proposition [L6] In particular, we see that the non-degeneracy 
condition combined with the Dirichlet boundary condition imposed on 7 yield equi- 
coercivity of the energies. The proof relies on an approach by G. Dal Maso et al. 
in [9] and combines the geometric rigidity estimate (see |12j ) with an estimate on 
7 that allows to eliminate free rotations (see [9] and Lemma [3~2l below) . Moreover, 
we prove a simple rigidity estimate for periodic variations (see Lemma l3.4|) . 

Section |4] is devoted to the proof of Theorem 11.11 and in Section [5] we present 
a linearization theorem that slightly extends results in [S] and is tailor-made for 
condition (W4). In Section [BJ we discuss the diagram (|1.3[) and prove Theorem 1 1.41 
Finally, in Section [7] we present two examples for which homogenization and lin- 
earization do not commute. 

Notation. We denote the space of real nxn-matrices by M™ and the subset of 
rotations, i.e. F € M n with F T F = Id and det F = 1, by SO(n). For matrices 
F,G6 M™ we define the inner product and induced norm by 

(F, G) := trF T G and \F\ := yJJF, F) respectively. 

Let k £ N. A measurable function u : R" — s- R is called fcF-periodic, if it satisfies 
u(y + z) = u(y) for almost every y 6 R™ and all z E kZ n . We define the function 
spaces 

L p pcr (kY) :={ue if oc (R") : u is fcF-periodic } , 

W^(kY) :={ue L? cl .(kY) : u E W^(R n ) } 

and likewise LP or (/cF;R n ) and W*g(kY;R n ). 

We denote the set of linear maps L from M" to M™ that satisfy (LA, B) = 
(LB, A) for all A,Be M" by T sym (rt), i.e. L € T sym (n) is a symmetric fourth 
order tensor. We associate to each quadratic integrand Q £ Q a unique map 
L e L~ r (F;T sym (n)) according to 

{H y)F, G) = Q(y^ + G)-Q(y,F)-Q( y ,G) ^ ^ p Q g ^ 

2. Stability of the non-degeneracy condition. Proof of Lemma [T73l 

In this section we prove Lemma IT~3l Let W € W(a,p). It is well-known that the 
homogenized integrand ■ M n — > [0, 00) is a continuous and quasiconvex map 

that satisfies the p-growth and p-coercivity condition (Wl) (see e.g. [TB], [S] and 
[3]). Condition (W2) is trivially fulfilled and Theorem O implies that has 
a Taylor expansion at Id, and therefore satisfies condition (W4). Thus, it remains 
to prove that the non-degeneracy condition (W3) is stable under homogenization. 
This is an immediate consequence of the following observation: 

Lemma 2.1. Let W : R™xM™ — > [0, +00] be a Caratheodory function Y -periodic 
in its first variable and suppose that W satisfies the non- degeneracy condition (W3). 
Then for all F G M" we have 

W^{F)>ldi^{F,SO{n)) 

where the positive constant a 1 only depends on the dimension n and the constant 
from condition (W3). 
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Proof. By definition we can find for each k € N a number € N and a map 
i>k e W^P{m k Y;R n ) such that 

wiZ } (F) + l>-^ I W(y,F + VMy))dy. 

k J 

We apply the non-degeneracy condition (W3) to the right hand side and obtain 
wiZ\F) + l>^ I d\st 2 (F + VMy),SO(n))dy. 

k J 

m 1" 

By replacing the map F \-> dist 2 (F, SO(n)) by its quasiconvexification Qdist 2 (-, SO(n)) 
we get a lower estimate: 

W£Z\F) + \ > \ [ ^ I Qdist 2 (^ + VMy),SO{n)) dy 



a \ m 



m k Y 



Because uikY is a quadratic domain and ipk is ra^K-periodic, we see that the inte- 
gral in the braces is bounded from below by Qdist 2 (_F, SO(n)) due to quasiconvexity. 
Thus, 

wj i Z\F) + l>lQdist 2 (F,SO(n)). 

K. Zhang proved in [53] that the quasiconvexification Qdist 2 (•, SO(n)) can be 
bounded from below by c„ dist (•, ( SO(n)) where c„ is a positive constant. Hence, 
we arrive at 

wiZ } (F) + l>^ist 2 (F,SO(n)) 
where a' := ^ > Q. Passing to the limit k — > oo completes the proof. □ 



3. EQUI-COERCIVITY BASED ON GEOMETRIC RIGIDITY 

In this section we prove Proposition 11.61 First, we would like to remark that 
Proposition ll.6l is an equi-coercivity statement. For instance the proposition implies 
that whenever a family of scaled displacements (g e ,h)s,h C L 2 (£l;M. n ) has equi- 
bounded energy, i.e. 

limsup I e ' h (g e , h ) < oo, 
( e ,/0-Ko,o) 

then it is relatively compact in L 2 (il;M. n ) and we can extract a subsequence that 
strongly converges to a map g £ _4 7 in L 2 (fl;M. n ). If additionally each element of 
the sequence {g £ ,h) has finite energy, then the construction of <F reveals that the 
relative compactness also holds with respect to the weak topology in W 1,2 (fl; K.™). 
In some sense this observation a priori justifies the presentation of the scaled energy 
in terms of the scaled displacement. 

A first step towards the proof of Proposition 11.61 is to show that the norm of a 
displacement gradient can be controlled by the associated energy. Because of the 
non-degeneracy condition (W3), we want to establish an estimate in the form 

(3.1) V 5 e A-y : J |V.g| 2 dx < C-^ J dist 2 (Id +h Vg(x), SO(n)) dx 

O f2 
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for a constant C that is independent of h. A key ingredient in the proof of this 
estimate is the geometric rigidity estimate by G. Friesecke, R.D. James and the 
first author: 

Theorem 3.1 (Geometric rigidity estimate [12 ). Let U be a bounded Lipschitz 
domain in K™, n > 2. There exists a constant C(U) with the following property: 
For each v £ W ' (17; R") there is an associated rotation R £ SO(n) such that 

J \Vv(x) - R\ 2 dx < C(U) J dist 2 (Vi;(a;), SO(n)) dx. 
u u 
Moreover, the constant C(U) is invariant under uniform scaling of U . 

In virtue of this result we can assign to each g £ A-y and positive parameter h a 
single rotation R £ SO(n) such that 

2 



ld+hVg(x) - R 



L 2 (0;I 



<C{n)^ J dist 2 (Id +hVg(x),SO(n)) dx. 



Let us assume for a moment that R = Id. Then the previous estimate and the 
non-degeneracy of W directly imply (|3.ip . In general a similar observation is valid: 
By taking the Dirichlet boundary condition imposed on 7 into account we can 
eliminate the free rotation. This has been shown by G. Dal Maso et al. in [5]. In 
particular, they proved the following result: 

Lemma 3.2 (see Lemma 3.3 in [9]). Let Q be an open and bounded Lipschitz 
domain in M™ and 7 a measurable subset of dQ with positive (n — 1) -dimensional 
Hausdorff measure. For F € M™ define 



\F\ 2 := min / \Fx - b\ 2 dH^Hx). 
1 'T fcGR" J 1 1 w 



Then there exists a positive constant C such that 



\F\ 2 <C\F\ 2 



for all matrices F that belong to the union of the cone generated by Id — SO(n) and 
the set of skew symmetric matrices in M™. 

In the lemma above T-L n ~ l denotes the (n — l)-dimensional Hausdorff measure. 
With this result at hand we are in position to prove the proposition: 

Proof of Provosition \l.b\ It follows from the proof of Proposition 3.4 in [9] that 

(3-2) < J d\ S t 2 (ld+h\7g(x),SO(n)) dx 

n 

for all g <E Ay and all h € (0, 1). Here and below c' denotes a positive constant 
that may vary from line to line, but can be chosen only depending on W and the 
geometry of 57 and 7. For the sake of completeness we briefly sketch the proof: 
In virtue of Theorem 13.11 we can assign to each g g A-y and h € (0, 1) a rotation 
R £ SO(n) such that 

J \ld+hVg{x) - R\ 2 < c' dist 2 (Id +h Vg(x), SO{n)) dx. 
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Since h Vg = (Ld+hVg - R) - (ld-R), we have 

2 



UsC'-WO < c '^2 / dist 2 (ld+hVg(x),SO(n)) dx + 2 



R-ld 



Set u{x) := (Id —R)x + hg{x) — uq where uq € R ra is chosen in such a way that u is 
mean value free. Because g{x) = on 7, we have u{x) = (Id— R)x — un on 7 and 
Lemma 13.21 implies that 



|Id-i?| 2 < C\ld-R\ 2 <C J \u{x)\ 2 dH^ix) < c' J \\7u(x)\ 2 dx, 

7 o 

where we used the continuity of the trace operator and Poincare-Wirtinger inequal- 
ity. Because of the identity V« = Id +h Vg — R, the right hand side is controlled 
by / n dist 2 (Id +hVg(x), SO(n)) dx and {EUl follows. 

By appealing to the non-degeneracy condition (W3) and Lemma |2. 11 we imme- 
diately see that 

mm{F' h (g),l£ om (g)} > c'-^ J dist 2 (Id +h Vg(x), SO(n)) dx 

o 

for all g € *4 7 and all s, h > 0. Since the energies are infinite whenever g ^ Aj, the 
previous estimate and (I3.2[) imply that 

imn{F' h (g),lb m (g)}>c'*(g). 

for all e,h > and g e L 2 (Q;K"). 

Next, we consider the energies (2S ) andX . It is not difficult to show that (W3) 
and (W4) imply that 

min I Q(y, F), Q^ m (F) } > d Isym^l 2 for all F G M™ and a.e. y € M". 

Thus, it is sufficient to prove that 

(3-3) Ibllw^cnjH") - c '/ l s y mV 3( x )| 2 dx for a11 9 G A- 

h 

This can be seen as follows: Let 5 g Ay. By Korn's inequality there exists a skew 
symmetric matrix if € M n such that 

llVff - K ||£ 2(n . M „) / |symV5(a;)| 2 dx. 


We set u(x) := g(x) — Kx — uq where tin € ffi™ is chosen in such a way that w has 
vanishing mean value. As before we can apply Lemma 13.41 and find 

\K\ 2 < C\K\ 2 < c' J \Vu\ 2 dx. 
Since V« = \7g — K, we find that 

IMIwi.*(n;R») - C ' ll V fflli 2 (0;M") < c ' ( II Vu IIl2 (0 ;M") + l#l 



Now (|3.3j) follows because the right hand side is controlled by ||sym Vg|| L2 (f2-M")- 

□ 
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Remark 3.3. In the proof we did not use the property that W satisfies the growth 
and coercivity condition (Wl). Moreover, we could deduce from Theorem 11.41 and 
general properties of T-convergence, that (lt ola ) and (If in ) and 1° are equi-coercive 
whenever (I £ ' h ) is equi-coercive. 

Another useful observation is the following: 

Lemma 3.4. There exists a positive constant c such that 

i f dist 2 {Id+hS7tp(x),SO(n)) dx > c J |W>| 2 dy 



kY 



kY 



for all h> 0, k e Z and maps ip G W*g(JtY; W 1 ). 

Proof of Lemma \3.4\ Set v{x):—x + hip(x). By Theorem 13. II there exists a rotation 
R G SO(n) satisfying 

2 , 



R-ld 



W>(y) 



dy<C 



dist 2 (Id+/iVV'(y), SO(n)) dy. 



kY kY 

The constant C is independent of h, k and tp, because the domain kY is obtained 
by uniformly scaling the cell Y. We expand the left hand side and see that 



i?-Id 



dy 



R-ld 



,'R-Id 
2( —r— , V^(y) 



|V^(y)r dy. 



kY 



k\' 



Because gradients of functions in Wp£(kY; R") are orthogonal to constant matrices 
(with respect to the standard inner product in L 2 (kY;M n )), we deduce that the 
integral over the coupling term in the middle vanishes and we immediately obtain 

|W(2/)| 2 dy < / dist 2 (Id +h^(y),SO(n)) dy. 



k-Y 



kY 



□ 



(4.1 
and 



4. Expansion of the multi-cell homogenization formula. Proof of 

Theorem 11.11 

For convenience we define for x S K™ and G £ M" \ {0} the remainders 

V^(a;,Id+G) - Q{x,G) 



r(x,G) 



rhom(G) :— 



II", 



(mc) 



horn 



(Id +G)-Q^JG) 



\G\' 



In order to prove Theorem ll.ll it is sufficient to show the following: For any sequence 
of matrices (G k ) in M" \ {0} with \G k \ -)• there holds 

(4.2) limsupr hom (G fe ) = 0. 

k— >oo 

Because the normalized sequence (\Gk\ 1 Gk) is relatively compact, it is sufficient 
to consider sequences that additionally satisfy 

Gk 



(4.3) 



H k := 



G k 



G 



as k — ^ co. 
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In the sequel we separately prove that 

(4-4) hmsup^V^Id+Gfc) < Q^JG), 

k->o \G k \ 

(4-5) liminf-l^^(Id+G fc ) > Q^G). 

Clearly, the validity of both estimates is equivalent to (|4.2j) . 

Step 1. We prove the upper bound estimate (|4.4I) . Because Q is a Caratheodory 
function quadratic in its second variable, the functional 

Wg(Y; K") 3 <p ^ Q(y, G + V<p(y)) dy 

is lower semicontinuous with respect to weak convergence and continuous with re- 
spect to strong convergence. Furthermore, the strict convexity of Q{x, •) on the 
subspace of symmetric matrices is sufficient to guarantee that the functional ad- 
mits a minimizer in W£g(Y;R n ) (see Remark [THJ). The inclusion C^ CI (Y;R n ) C 
Wp£(Y;M. n ) is dense; thus, by the strong continuity of the functional we find for 
every 77 > a map i\i £ C^ r (Y; E n ) such that 

= mm f Q(y,G + V<p(y)) dy > [ Q(y,G + W^(y)) dy - v . 

Y Y 

Based upon this choice we derive an upper bound for the left hand side in (|5.5p : 
By construction the multi-cell homogenization Wu^f (•) is bounded from above by 
the one-cell homogenization w£2n(")i thus, we obtain 

w£^(Id+G fc ) < Will(ld+G k ) < J W{y,ld+G k + \G k \V^(y)) dy. 

Y 

Inequality (*) follows from directly follows from the definition of W^L. We expand 
the integrand on the right hand side and deduce that for almost every y £ Y we 
have 

-±^W(y,ld+G k + \G k \ V^(y)) < Q(y, H k + W(y)) + r(y, G k + \G k \Vi>(y)) 
\Gk\ 

where the remainder r is defined according to (14.11) . Because (Gk) vanishes and 
t/> £ W^iYiW 1 ), we have 

G k + \G k \V^(y)^0 uniformly; 

thus, condition (W4) implies that 

limsup / r(y, G k + \G k \ Vip(y)) dy = 0. 

k— >oo J 
Y 

Consequently, the previous estimates, the convergence H k — > G and the continuity 
of Q lead to 

l imsup _L^)( Id+Gfe) < f Qfa G + v ^(y)) dy < Q$ m (G) + T). 

fc— >oo \(-Jk\ •* 

Y 

Because 77 > can be chosen arbitrarily small, (|5.5j) follows. 
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Step 2. We prove the lower bound estimate (|4.5[) . We only have to consider 



the case where 



1 



k^oo | Gfc | 

is finite. We pass to a subsequence, that we do not relabel, such that 



lim inf ■ 



1 



,2 "hom 



1 



(Id +G k ) = lim sup — l,^' (Id +G fc ). 



fc ^o | Gfc |" — ' |G fc | 
By definition, for all k 6 N there exist a number nik 6 N and a map 

such that 



with 



(4.6) ^)(Id+G fc ) + ^>^ 



W(y,U+G k + \G k \VMy)) dy. 



This suggests to establish the lower bound estimate by applying the expansion in 
condition (W4) to the integral on the right hand side. Clearly, if \Vip k (y)\ was 
bounded by a constant independent of y and k, condition (W4) would immediately 
imply (|4.5j) . However, the sequence (V^fc) is only bounded in the following sense: 

I 



(4.7) 



G := lim sup ■ 

k— >oo m i 



\Vipk\ dy<oo, 



as can be seen by appealing to the non-degeneracy condition (W3) and Lemma[ 
Therefore, we distinguish points y £ rrikY where \^tpk{y)\ is sufficiently small from 
those where |VV>/c(y)| is too large for an expansion. More precisely, we define the 
set 



(4.8) 



Y k :={ye m k Y : |W> fc (y)| < |G fc p 1/2 } 



and denote the associated indicator function by Xk- The proof of the lower bound 
estimate is divided in two steps. First, we show that 



(4.9) 



lim inf i w££ (Id +G k ) > lim inf -L f Q (y, Xk (H k + V Vfc)) d y 



and in a second step, we prove that 
I 



(4.10) 



lim inf ■ 

k-^oo rri'i 



Q{y, Xk(H k + W*)) dy - Q^JG) > 



It is obvious that the combination of both estimates justifies (|4.5p . 
Step 3. (Proof of (HJ])). By construction we have 



ess sup 

yen 



G k + \Gk\VMv) 



and in view of condition (W4) we see that 
(4.11) 



lim sup 

fc^o ml 



Xk 



W(y,Id+G k + \G k \W^ k ) 



\G k \ 







Q(y,H k + V^ k ) 



dy = 0. 
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Now the non-negativity of W implies that W(y,F) > Xk{y)W(y, F) for almost 
every y £ irikY and all F € M"; thus, estimate (|4. 6[) and ()4.11j) immediately imply 
that 



k + V?/>fc) dy. 



"^^(Id+G^^liminf-L / XkQ{y ,H, 
Because \k takes only values in {0, 1}, we see that 

Xk(y)Q(y,H k + VMv)) = Q(y, Xk{y){H k + VMv)) 

for all y £ m k Y and (|4.9[) follows. 

Step 4- (Proof of (|4.10jl ). This is the heart of the matter. Obviously, if the 
indicator function \k was equal to 1, then the integral in (|4.10[) would be bounded 
from below by Q?' (ffs) and the estimate would follow from the continuity of 
Qhom- The g enera l case would follow if we knew that 



l - [ (1 - X k)\H k + W^ k \ 2 dx^O. 



l k Jm k Y 

Because a priori IVf/'fe) 2 could concentrate on the set where \ k — 0, this is not 
obvious at all. Since we aim for a lower bound, we expand Q{y, H k + V^fc) around 
the minimizer G + V"0g (see below) . Then the most dangerous quadratic term has 
a sign. 

Thus, let Vg € W^(Y;R n ) satisfy 



Q$JG)= / Q(v,G + VMy)) dy- 



We extend ipG by periodicity to 1™. Since Q is F-periodic and convex, it is not 
difficult to check that ipc is a ls° a minimizer of the functional 

(4.12) W$(m k Y',R n ) B 1> » J Q{y,G + V^(y)) dy. 

Because Q(y, •) is a quadratic form, the inequality 

Q(y, A) - Q(y, B) > 2 (h(y)(A - B), B) 
is valid for all A,B£ M" and almost every y. We apply this inequality with 

A = Xk(y)(H k + VMv)) and B = G + Vip G (y). 
Now integration over mkY leads to 
1 



2m£ 



Q(y, Xk{H k + Vipk)) - Q{y, G + Vi/>g) dy 



m k Y 



>\ ( (h(y)\(xk(H k + V^k)-(G + V^ G ) 
ml J \ L 



, G + V^g) dy. 



m k Y 
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The right hand side can be rewritten as 

\ f (h(y)(H k -G), G+Wg) dy 
ml J 

m k Y 

+\ [ {L(»)(V^fc-VVc),G+V^ G )d» 
m n k J 

~ f (L(j/)(l-»)(ft+V^),G + Vfe)dy 

m k J 



m k Y 

.,-(1) , r(2) , r(3) 
' k ~t~ 1 k ^ 1 k ■ 



In the following we prove that all three integrals vanish as k — ¥ oo. We start with 

l k 



the first integral j£ . Due to the ^-periodicity of L and ipQ, we have 



jf } = / (L(y)(F fc -G), G+V^g) dy 



and (|4.3I) implies that — > as fc — >• oo. We consider the second integral 
If,'. It is easy to check that Ii ' is exactly the Euler Lagrange equation of the 
quadratic minimization problem associated to the functional (|4.12[) . Because tpa is 

(2) 

a minimizer and the map y ip k — WG an admissible test function, we have r k ' = 
for all fc G N. 

(3) 

Finally, we consider the third integral If, . By applying the Cauchy-Schwarz- 
and Holder-inequality, we find that 
(4.13) 

2 



r(3) 
1 k 



\ m k Y / \ 



" ''\— I |///, + Vr A | J <ly 1 | ^ J \(l- X k)(G + ^ G )\ 2 dy 

m k Y 



where c' is a positive constant independent of k. In virtue of (|4. T[) . it is sufficient to 
prove that the second integral vanishes as k — > oo. By construction the multi-cell 
mkY is the disjoint union of the mjj translated cells Y +£ with ( e := mkYC\1 n . 
Consequently, we can rewrite the second integral in (|4.13l) according to 



-4 / l(l-Xfe)(G + V^ G )| 2 dy 
ml J 



"k 

m k Y 



dy 



where we used the fact that the map y h-> (1 — Xfc) only takes values in {0,1}. 
Because the map i\)q is K-periodic, the right hand side simplifies to 

[ Xk(y)\G + V^ G (y)\ 2 dy with Xfcfa) := (l - Xk(» + £))• 
v mfe «ez fc 
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Recall that %fc denotes the indicator function of the set Y k defined in (|4.8|) . By 
definition there holds f < \Gf-\ |V?/>fc(2/)| for a.e. y € m k Y\ Yk. Thus, we estimate 

f\Xk(y)\ dy< / (1-Xk(y))dy<-L [ \G k \ |V^(y)| 2 dy < C \G k \ 
J m k J ml J 



m k Y m k Y 



where C denotes the constant from (|4.7p . Consequently, the sequence Xk strongly 
converges to in We claim that 

(4.14) X k \G + V^ G | 2 -> strongly in L\Y). 

Since Xk ^ in L 1 (y), the left hand side converges to in measure. On the 
other hand, by construction we have < < 1, and therefore the left hand side is 
dominated by the map \G + Vipc(')\ 2 which belongs to L 1 (l r ). Thus, by dominated 
convergence assertion (|4.14p follows and ijp vanishes as k — > 0. So far we have 
shown that 

\ [ Q{y,Xk{H k + V^ k ))-Q{y,G + V^ G {y))dy^Q 
ml J 

m k Y 

as k — > oo. In virtue of the Y-periodicity and convexity of Q, we have 
QhL(G) = ~^ I Q(y,G + ^ G (y))dy 

" l k J 

m k Y 

and (|4~10| follows. □ 

5. Linearization. 

G. Dal Maso et al. proved in [5] that linearized elasticity can be obtained as 
a r-limit from nonlinear, three-dimensional elasticity. The following theorem is a 
variant of their result adapted to assumption (W4) . 

Theorem 5.1. Let W : fixM™ — > [0, oo) be a Caratheodory function and suppose 
that 

(KU v \W(x,ld+G)-Q(x,G)\ 
(5.1) hmsup —2 = 0, 

G^O, \G\ 

where Q : J!xM" —¥ [0,oo) is a Caratheodory function quadratic in its second 
variable and bounded in the sense that 

VG € M" : ess sup Q(x, G) < c\ |G| 2 
for a suitable constant c\ > 0. We consider the functional 

W(x,Id+hVg(x)) dx if g € Ay 



£ h (g) 



n 

oo else, 



and assume that there exists a positive constant ci such that 

(5.2) £ h (g) > c 2 ||ff||^i, a(n . R „j for all h > and g G A, 
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Then the family (Eh) T-converges with respect to strong convergence in L 2 (il;R") 
to the functional 

— f Q(x,Vg(x)) if g G Ay 

Q 

oo else. 



£lin(ff) 



h 2 



Remark 5.2. Condition (|5.2|l implies equi-coercivity of the functional (£ h ) in 
L 2 (Q;R"). In virtue of Proposition 11.61 we see that the combination of the non- 
degeneracy of W and the Dirichlet boundary condition is a sufficient condition for 

in. 

Proof. Since the strong topology of L 2 (Cl; R") is metrizable, we can use the sequen- 
tial characterization of L-convergence. Thus, we have to prove the following: 

(1) (lower bound) for every g G L 2 (f2;R n ) and every sequence (gh) converging 
to g in L 2 (fl:R n ) it is 

liminf £ h (g h ) > £u a (g). 

h—tO 

(2) (recovery sequence) for every g G L 2 (il;M. n ) there exists a sequence (<??,,) 
converging to g in L 2 (f2;R") such that 

lim £ h (g h ) = £iin(g)- 

h— >-0 

Step 1. (Recovery sequence). We only have to consider the case g G A 7 . By 
assumption the inclusion ( W 1,0Q (f2; R") n Ay ) C Ay is dense; thus, there exists a 
sequence (gh) in W^ 1,oc (57; R") n Ay that strongly converges to g and that satisfies 

(5.3) ess sup | VflVi (a;) | < — ^= for all h > 0. 

Thus, by (jSTTj) we get 
1 

h 2 



(5.4) / W(x,U+hVg h (x)) dx = I Q(y, Vg h (x)) + v(x,h Vg h (x)) dx 



with 



lim 

h-yO 



v(x, hWgh(x)) dx 



= 0. 



Since the quadratic integral functional in (|5.4[) is continuous with respect to strong 
convergence and gh — > g strongly in W 1,2 (fl;W), we see that 

lim £ h (g h ) = £v m (g)- 

Step 2. (Lower bound). Let g G L 2 (Q;W l ) and (gh) a sequence converging to g 
in L 2 (n;R n ). We show that 

(5.5) liminf £ h (g h ) >£i in (g). 

h— >0 

As usual it is sufficient to consider the case where 

liminf £ h (gh) = limsup £ h (gh) < +oo. 

In this case we can assume without loss of generality that each gh belongs to the set 
Ay. The equi-coercivity condition on (£ h ) (see (|5.2p ) implies that (gh) is relatively 
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compact with respect to weak convergence in W 1,2 (Q; W 1 ). As a consequence we 
can pass to a subsequence (that we do not relabel) such that 

g h ^g weakly in W 1)2 (fi; K n ). 

Because A 7 is a (seq.) weakly closed subset of W 1,2 (f2;R n ), we find that g e Ay. 



i 7 

We define the set 



Q h :={xen : \Vg h {x)\ < h^ 2 } 



and denote the corresponding indicator function by Xh- Since (Vgh) is a bounded 
sequence in L 2 (f2; M"), we can estimate the measure of fi \ Q/j according to 

(5.6) \n\Q h \<h J \Vg h {x)\ 2 dx <hC 

for a suitable constant C . Now we establish a lower estimate by utilizing the 
expansion of W for points x € flh and the non-negativity of W for points x € Q\Qh- 

W(x,ld+hVg h (x)) dx > / Q(x, Vg h {x)) + t(x, h Vg h {x)) d x. 



bmce \hVg h {x)\ < Vh for all x G Vih the properties of the quadratic expansion of 
W lead to 

liminf / W(x,Id+h Vg^fx)) dx > liminf / Q(x, Xh(x) Vgh(x)) dx. 

n n 

Since V.g^— ^Vg weakly in L 2 (f2;M") and x/i — ^ 1 boundedly in measure (due to 
dSHO, we find that 

Xh Vg h ^ Vg weakly in L 2 (ft; M") 

and the sequentially weak lower semi-continuity of the quadratic integral functional 
leads to the desired lower bound estimate (15.51). □ 



6. HOMOGENIZATION AND LINEARIZATION COMMUTE 



In this section we prove that the diagram (jTT3J) in Theorem 11.41 commutes. We 
have to show that the T-convergence statements 

(1) P^Ii, (2) T- h AiL 

(3) ig n -A i° (4) t ra Ar° 

hold w.r.t. strong convergence in L 2 (f2; M"). 

First, we discuss the r-limits corresponding to linearization, i.e. (1) and (4). To 
this end we recall that Proposition 11.61 implies that 

(6-1) mm{T> h (g),l£ om (g)}>c\\g\\ 

for all g E A~ t and h,e > 0. As a consequence we see that convergence (1) directly 
follows by applying Theorem 15. II to the sequence (I e > h )h- 

For the justification of convergence (4), we first apply Theorem ll.il and see that 
admits a quadratic Taylor expansion. It is well-known that the homoge- 
nization Q^ m of a quadratic integrand Q € Q is a quadratic form over R", and 
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therefore satisfies the condition of Theorem 15.11 Moreover, Proposition II .61 proves 
that the sequence (Z^ om ) is equi-coercive and (4) follows by applying Theorem 15. II 
The T-limits corresponding to homogenization, i.e. (2) and (3), follow by stan- 
dard results. More precisely, convergence (3) is a classical result in convex homog- 
enization (cf. e.g. [10], [18] and note that these results can be restated in the 
language of T-convergence). The non-convex homogenization, convergence (2), can 
be justified by the following result: 

Theorem 6.1 (A. Braides [4], S. Miiller [16]). Let f2 be a Lipschitz domain in R". 
Suppose that W : R™xM™ — > R is a Caratheodory function that is Y -periodic in 
its first variable and that satisfies the growth and coercivity condition (Wl) with 
p G (1, oo). We consider the Junctionals 



£*(u) :-- 

and 

£hom(u) 



o 



W(x/e, Vu(a:)) da; if u € W hp (n; R") 

oo else, 

W^(Vu(x)) dx ifu€ W 1 ' p {il;R n ) 

else. 




Then (£ 6 ) T-converges to £hom with respect to strong convergence in L p (Q:R n ) as 



The previous theorem can be adapted to the boundary conditions considered in 
diagram (|1.3[) . More precisely, one can show that 

for every u € W 1,p (tt; R n ) there is a sequence (u s ) such that u e — >• u 
strongly in L p (il;R n ), u e — u on dil and lim£ e (u e ) = £t om (u). 

e— >0 

This can be seen by a standard gluing argument (cf. e.g. [16[ IS"]). In particular, 
we can apply this modification in the case p > 2 and to limiting deformations 
u £ M /1,2 (0;R n ) satisfying u(x) = x on 7. Hence, the theorem implies that (I e - h ) £ 
T-converges to Ijjom w ith respect to the strong topology in L 2 (f2;R"). 

In conclusion we see that both paths (1), (3) and (2), (4) lead to the same 
limiting functional 1° and the diagram commutes. 



7. Failure of commutativity 

In this section we present two examples illustrating that homogenization and 
linearization cannot be interchanged in general. In the first example we consider 
energy densities of class W(a,p) and show that the commutativity might fail if the 
expansion is not centered at identity. Secondly, we elaborate on the importance of 
condition (W2) and (W3), which imply that the material has a single energy well 
at SO(n) with quadratic growth. As an example we discuss a perforated composite 
with a prestressed component and give an indirect argument that homogenization 
and linearization do not commute at identity — although the homogenized energy 
density is stress free at identity and satisfies Wj"^(Id+G) = 0(|G| 2 ). 
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7.1. Counterexample I: Non-commutativity for expansions away from 

SO{n). In this section we argue that in general the commutability of homoge- 
nization and linearization does not hold for expansions centered at F ^ SO(n). 
Roughly speaking, the reason is the following: If F ^ SO(n) we may find a nonlin- 
ear material W G W(a,p) with wj^(F) < Wy^ m (F), which means that long- wave 
oscillations with period ke, k > 2, lead to limiting energies that are lower than those 
obtained by "one-cell" oscillations with period e. On the other side, the linearized 
energy and its homogenization is always quadratic, and due to the non-degeneracy 
of the material stable in a neighborhood of SO(n). Thus, low energy states related 
to long- wave oscillations are ignored by the homogenization of the linearized energy; 
and therefore homogenization and linearization do not commute for F SO(n) in 
general. 

In the rest of this section we present an implicit, but rigorous formulation of 
the idea above by varying an example introduced by the first author in [IB]. Let 
Wo : M 2 — > [0, oo) be a frame indifferent integrand of class C 3 . We suppose that 
Wo satisfies the local Lipschitz- and growth condition (Wl) with p = 2, and 

(7.1) - dist 2 (F, 50(2)) < W (F) < cdist 2 (F, 50(2)) for all F G M 2 . 

c 

We consider the periodic stored energy function 

W(y, F) := ( X (y) + a(l - X (v)) )W {F), 

where a is a small positive parameter and \ denotes the indicator function of the 
periodic pattern 

P := { (£1,2:2) £l 2 : 3k G Z such that x 2 E [k, k + |) | 

Thus, W describes a composite material with a layered microstructure consisting 
of a stiff component on P and a soft matrix with stiffness a. 

As a main result we prove a necessary condition for the property that lineariza- 
tion and homogenization commute for all expansions in an neighborhood of Id: 

Proposition 7.1. Let U be an open neighborhood of Id. Suppose that for each 
expansion at F e U homogenization and linearization commute in the following 
sense: There exists fff € M 2 such that 

VG G M 2 : W^(F + G) = W^ 3 (F) + (a F , G) + Q[ om (G) + o(\G\ 2 ) 

where Ohom denotes the homogenization of the guadratic integrand 

Q F {y,G) :=±D 2 W(y,F)(G,G). 

Then 

a > c(W Q ,U) 

for a positive constant c(Wq, U) > that only depends on Wo and the neighborhood 
U. 

Remark 7.2. Theorem 11.11 implies that the expansion holds for F — Id with up = 0. 

Remark 7.3. The result can be read as follows: For any open neighborhood U 
of Id we can find a material (e.g. by choosing a sufficiently small) such that 
homogenization and linearization do not commute for all expansions centered in U. 
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The result follows from a slightly stronger statement that we prove below. In 
particular, it is sufficient to study the response of the homogenized material to 
compressions of the form 

x h> F$x, F$ := Id-5(ei(g)ei), 6 > 

For future reference we set 

Q 6 (y,G) :=^D 2 W(y,F s )(G,G) 

and let Q^ om denote its homogenization. 

In Q15] it was shown that if 5 > and a is sufficiently small, the stiff part of the 
material, which resembles an ensemble of aligned rods, starts to buckle and allows 
to bound the limiting energy linearly in a. In contrast to this, we prove that the 
non-degeneracy of W implies that Qt om is stable for small S. Both observations 
can be quantified as follows: 

Lemma 7.4. There exist positive constants Sq and cq (depending only on Wq) with 
the following properties: For all < S < Sq and a > we have 

(7-2) QhomCeiOer) > -, 

CO 

(7-3) W^Z\Fs)<c a. 

We postpone the proof to the end of this section and present a quite immediate 
consequence of the lemma that already implies the validity of Proposition 17.11 

Lemma 7.5. Define the map 

/:R->R, f(S):=w£$(F B ). 

Then for almost every S € [0, 5q] the map f admits a quadratic Taylor expansion of 
the form 

(7.4) f(5 + A) = f(S) + a s \ + q s X 2 + o(A 2 ) 

for suitable numbers as G R and qs > 0. 

(1) If 5 = 0, the expansion (|7.4[) is valid for 

(T = and q a = Q° m( e l® e l)> 

i.e. linearization and homogenization commute at Fq = Id. 

(2) Suppose that linearization and homogenization commute at Fg for all a.e. 
S e [0,6 a ], i.e. 

QS = <3hom( e l® e l) f or a - e - $ G [0,6o]. 

Then 

a > 

Remark 7.6. In the lemma above the constants <5o and cq only depend on Wo. 
Thus, if the soft matrix of the composite material is sufficiently weak, i.e. a « 1, 
statement (2) suggests that homogenization and linearization may only commute 
in a very small neighborhood of Id. 
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Proof of Lemma |7.5[ By Lemma 11.31 the map W^™^' : M 2 — > K is continuous and 
quasiconvex, and therefore rank-one convex. Because Fg — F\ is always a rank-one 
matrix, we deduce that / is convex and continuous. Thus, Aleksandrov's Theorem 
implies that / admits an expansion of the form (|7.4[) for a.e. 5 € [0,<5o]. Statement 
(1) is a direct consequence of Theorem 11.11 applied to G = — A^tgiei). 
We prove (2). By assumption we have 

/"(*) = Qhom^er) 

for a.e. S. We apply (|7.2I) from Lemma \7. 41 and find that /" > Cq 1 almost every- 
where. Thus, the map 

9(5) ■= f(S) - ^-S 2 

2c 

is convex and by part (1) we have <?'(0) = f'(0) = 0. Thus g attains its minimum 
at with g(0) = 0. By QfJb we S et 

7^-S 2 < f(S) < c a for all S € [0, 6 ] 
2c 

and the estimate for a follows. □ 



Proof of Lemma \U\ In the sequel c', c" denote positive constants that may change 
from line to line, but can be chosen only depending on 8$ and Wo- 



Step 1. We prove (|7.2p . Let Qq denote the quadratic term in the expansion of 



Wo at Fs € M 2 , i.e. 



Q S (G) :=±D 2 W (F S )(G,G). 



Because F = Id and Wg is frame-indifferent and non-degenerate, the estimate 

Q°o(G)> ^|symG| 2 

c 

holds for all G € M 2 . The map Wo is of class C 3 , and therefore a standard 
perturbation argument shows that 

Q S o(G) > - |sym G| 2 - c"6 |skw G| 2 
c 

for all 5 £ [0,<5o] where So denotes a small positive constant that only depends on 
W . 

Let Yo := [0, l)x [0, ^) denote the reference domain of the stiff component and 
consider the function space 



V := {i> e W^iYo;^ 2 ) : ^{0, x 2 ) = ^(1, x 2 ), / ^)dx = o}. 

Note that for all tp € V a Korn inequality of the form 

/ IskwVf/'l 2 dy < c / IsymVV'l 2 dy 
Jy Jy 

holds. As a consequence we deduce that 

m{8o) := inf inf / Q 5 a{ex®ex + Vt()(x)) dy > 
<5e(o,<5 )iAev Jy 

provided So is sufficiently small. 
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Next, we are going to show that the estimate above yields a lower bound for 

•l, 

pc 



Q 5 hom (G): Let 6 G [0,6 ] and $ G W^ 2 (Y;R 2 ). Then 



(y, (ei<8ei) + V^) dy > / Qg((ei®e x ) + W) dy > m(S ) > 

and (|7.2I) is valid for all Co > m(5o) _1 . 

Step 2. We prove (|7.3p by constructing a sequence (ufc^gN such that 

(7.5) liminf / W (kx, Vzifc(x)) dx < c a 

and u k ^F s x + ip(x) weakly in W 1 ^ 2 (Y;R 2 ) with tp £ W^ 2 .(Y;R 2 ). Because the 
functional 

W h2 (Y;R 2 ) 3 u i-> J W (kx, Vu(x)) dx 
T-converges in the strong topology of L 2 (Y;M. 2 ) to 

W 1 ' 2 {Y;R 2 ) 3 u ^ J wfi(Vu(4)) dx 
as k — > oo, (|7.5p implies that 

(7.6) / W^(F S + V<p(x)) dx < liminf / W (kx, Vu fe (x)) dx < c a. 

The homogenized integrand 11^™' is quasiconvex, and therefore the left hand side 
is bounded from below by W^™^(Fs). 

It remains to construct the sequence (uk). We follow the idea in [16]: The 
microstructure of the material can be regarded as an ensemble of thin rods, aligned 
to the ei-direction and embedded in a soft matrix. The subsequent construction 
realizes the compression x i— > F$x (more precisely, a periodic variation) by bending 
each rod. In this way, locally the deformation in the stiff component is close to a 
rigid motion, while large strains only evolve in the soft matrix. 

For the precise construction, let vg G C°°([0, 1]; R 2 ) denote a curve, parametrized 
by arc-length and rendering an arc with length 1 connecting the points (0, 0) and 
(1- (5,0), i.e. 

v s (0) = (0,0), « tf (l) = (1-5,0), \v' s (x 1 )\ = l, K(xi)| = /t 4 

for some k$ > 0. Let ng G C°°([0, 1]; R 2 ) denote a normal field associated to the 
curve, so that the map x\ <-> R$(xi) := v' s (xi)<g>ei + ns(x\)®C2 only takes values 
in 50(2). 

For x € Y define 

' x 2 e 2 +v 5 (x 1 ) + (x 2 - |) (n s (x 1 ) - e 2 ) if x 2 G [|, 

for some i G No, 



Wfc(x) 



x 2 e 2 + -x 2 )(n (5 (x 1 ) -e 2 ) if ^ G [2|tl, i±l) 

for some i £ No. 



By construction we have Uk G W 12 (Y; R 2 ) and Uk — > vs + x 2 e 2 uniformly. The 
sequence (uk) is bounded in W 1 ' 2 (Y; R 2 ) and therefore 

Uk^F s x + tp(x) weakly in T4^ 2 (y; R 2 ) 
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where ip(x) := vg(xi) + xiei — Fgx. Furthermore, it is easy to check that (p can be 
identified with a map in Wp£(Y;R 2 ). 

Vw„ is close to a rotation on the strong part of the material. More precisely, we 
have 

,.2 * 



dist a (V« fc (a;), 50(2)) < ( '-j- ) k~ 2 for all x 2 e [£, ^) 



while on the soft part the deformation is at least bounded, i.e. 

\Vu k {x)\ 2 <c' for all x 2 e 
Thus, by using the definition of W and the growth condition (|7.1[) . we see that 

fc-l / fl /.(2i+l)/2fc 



W(kx,Vu k (x)) dx < cj^l / / dist 2 (Vw fe ,50(2)) da; 2 

+ a / / |Vu fe | +2 dx 2 dxi 

JO J(2i+l)/2fe / 

< c'(fc- 2 + a) 

and (17.61) follows as k — >• oo. □ 



7.2. Counterexample II: Non-commutativity at identity for a perforated, 
prestressed composite material. In this section we construct an example in 
dimension 3 for which the commutativity of linearization and homogenization at 
identity fails, although the identity is a natural state of the homogenized material. 
The material under consideration is a composite with a prestressed component 
that violates property (W2). Furthermore, the material is perforated and the non- 
degeneracy condition (W3) is violated in cc3-direction. 

A commutativity statement in the form of Theorem ll.il comparing the homoge- 
nized and the spatially heterogeneous energy density, is not suited for a prestressed 
material due to the lack of a stress free reference configuration. Therefore, we 
introduce a weaker property that is motivated by the following observation: 

Corollary 7.7. Let W € W(a,p). ForkeN let 



(F) :=inf J i_ J W(y,F + V<p(y)) dy : <p G Wg{kY;R n ) 



kY 

denote the k-cell homogenization ofW. Then 

(7-7) lim^^dd+Z.G) = Jta^eOd+ftG). 

for allkeN and all G € M n . 

Proof. The statement is implicit in the proof of Theorem ll.il Indeed, in Step 1 we 
showed that limsup 

W W hllu( ld + hG ) < QhoL( G )> while in Ste P 2 we Proved the 

inequality lim inf fcW™ (W +hG) > Q^JG). Since < <^ < w£ m we 

can draw the conclusion. □ 
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The corollary indicates that in situations when commutativity holds, close to 
identity relaxation over a single cell is sufficient for homogenization. Proposition l7.ll 
shows that if we replace Id by a matrix F $ SO(n), identity (|7.7[> fails for some 
materials of class W(a,p). This suggests to introduce the following property: We 
say that homogenization and linearization in direction G € M" commute at Id, if 

for all k G N there exist cr (fc) , a € M n such that 
(Co) ^( W ^+hG) - W^ (Id) -A(a«, G )) 

= it h ( w ™ (Id +hG) ~ ^™ (Id) ~ h g> ) 

Note that we distinguish between the directions +G and —G. By Corollary 17.71 
materials of class W{a,p) satisfy QC G D for all G € M". 

In the following, we present an example in dimension n — 3 that violates QCg| ). 
Let us first describe the geometry and the energy density of the composite. The 
composite's geometry is given by two subsets of the reference cell Y := [0, l) 3 : 

Y := [0,l) 2 x[0,i) and Y p := [0, l)xB p , 

where B p := {(y 2 ,y 3 ) ■ (ite - + (ife - f ) 2 < p 2 } , < p « 1. 

The set Y p is a cylinder with center line {(j/i, \ , f ) : J/i € [0, 1) } and small radius 
p. Let Wo be of class W(a, 2). We suppose that Wo is frame-indifferent, i.e. 
W {RF) = W {F) for all matrices F G M 3 and rotations i? G SO(3). Set 5 := 
(Id+s(ei®ei))- 1 with < s < |. We denote by W : R 3 xM 3 ->■ [0,oo) the 
K-periodic energy density 

(W Q (F) ifyeFo, 
W(y,F):= }\Vo(FS) if y G Yp, 

[o if yer\(r ur p ). 

Clearly, for the material occupying Y"o every rotation is a natural state, while on Y p 
the material is stress free only for matrices F = RS~ X with R G S'0(3). 

Lemma 7.8. The following properties hold: 

a) W is frame-indifferent, has quadratic growth at infinity and is non- degenerate 
in X\- and X2-direction in the sense that 



(W3') 



J W{y, F + Vif>(y)) dy > c dist 2 ^, SO (3)) 

for all 4) G W p ^ 2 (r;M 3 ) with [ [ d 3 tp d Vl dy 2 = 0. 



o Jo 



b) For all k eN we have inf FGM 3 W^(F) > 0. 

c) (Natural state). 

W£?(H) = 0= Egn W^(F). 

d) (Expansion at Id). For all G £ M 3 we ftave 



limsup— W^Id+W?) < oo. 
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e) (Failure of \Cg\ )- There exists G £ M 3 such that for all k £ N and all 
cr(fc) e M 3 with 



(7.8) 



lim sup -pr 



we have 



W^ m (Id+hG)-W^ a (Id)-h/a^, G 



< oo 



(7 - 9) ^ h (WwCM+W?) - <L(Id) - h G 



1 



In a nutshell, the idea behind the construction is the following: Close to identity 
and for finite k, the prestressed component contributes a positive amount of energy 
and yields a shift of the natural state to some F £" SO(3). When k increases, 
the aspect ratio (thickness / length) of the cylindrical component and (as a conse- 
quence) the contributed elastic energy decreases. We show that in the limit k — > oo 
the energy contribution of the Yp-component vanishes. Furthermore, we argue that 
vanishes for short maps. 

Proof of Lemma \7.8\ Let us introduce the functionals 



W(y,F- 



kY 



1 

i 

k 



W (F- 



W, 



(('■ 



-Vp)) dy, 
Vip) dy, 

■V(p)S) dy, 



(0,fc) 2 x(0,±), 



(0,k)xB p . 



and the function spaces 



X 



X, 



(fe) := 
(*) ._ 



{ 



pe 

I Zjp. 



<p £ I« } 
</> € X< fc > } 



inf 



vex 



Note that by definition we have W^^(F) 

The proof is divided in several steps. In Step 1 we argue that since the com- 
ponents Y"o and Y p are separated, we can split 1^ into 1^ and 1^ . In Step 2 
we argue that Ip k \-', F) is positive for F close to Id and show that its contribution 

(k) 

vanishes as k — > oo. In Step 3 we analyze the behavior of Iq . In particular, for the 
asymptotic behavior k — > oo we appeal to a T-convergence result for membranes 
(see [21[7]). In Step 4 we draw the conclusion. 

Step L (Splitting of /W). We claim that for all (<p, ip) £ ifxlf there exists 
$ £ XW such that 

(7.10) /« ($; F) = I { k) (if, F) + jW F) 

and 



(7.11) 



inf / (fe) ($;F)= inf (cp; F) + inf I^(ip;F) 
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Statement (|7.10[) easily follows from a periodicity property of the composite's mi- 
crostructure. Namely, we can rewrite the subset of the multi-cell kY occupied by 
the Yo-material as a union of k translations of the set Zq ; and the subset occupied 
by the Yp-material as k 2 translations of the set zf^ : 

y^--= U (t+ Y o)= U ((o,o ;l ) + 4 fe) ), 

Cez3n[o,fe) 3 iezn[o,fc) 
^ez3n[o,fc) 3 ijezn[o,fe) 



For y G F (fe) U Y p (fe) set 



' \p{y) if y G (0, 0, i) + Z^ fe) for some i G Z 
•0(y) h y e (0, + Z^- 1 for some i, j G Z 



Since Y^ n Yp fc ^ = we can extend $ to a fcY-periodic map $ G X' fc ). Since 
W(y,-) = whenever y G (fcY") \ (y" (fe) U Y p (fe) ) identity (T7TTU1) follows. (17111) 
follows from (|7.10p and the observation that ">" trivially holds in (|7.1ip . since the 
minimization problems on the r. h. s. are less constrained. 
Step 2. (Estimates for I { p k) ). Set 

We claim that 

(7.12) Vfe G N, G G M 3 : liminf m( fc) (Id +/iG) > 



(7.13) VG G M 3 : lim sup m( fe) (Id +/iG) - , , , 



h->0 

'O ifGe x = 



Proof of (|7.12p by contradiction: Suppose there exists a vanishing sequence of 
positive numbers (hj)j^jq and a sequence (ipj) C Xp fe ^ such that 

(7.14) lim inf (tpj ; Id +h s G) = 0. 

Consider the deformation Uj(x) := (Id +hjG)x + il>j{x) and the mapping Wj(z) := 
Uj(Sz) defined for z G S~ l Zp . Since Vtty(z) = Vuj(Sz)S, a change of variables 
yields 

(TT^fc / s _^ d2 - 1 4> Wo(V«,(»)5) d* = /(^. ; i d+/l .G). 

Thus, by (|7.14l) we have (up to a subsequence) 

(7.15) —}— [ W o (V Vj )dz^0 (asj^oc). 

In combination with the non-degeneracy condition (W3) and by appealing to geo- 
metric rigidity (see Theorem 13. lj) we find (after passing to a further subsequence) 
that there exists a rotation R G 50(3) such that 

(7.16) Vv 3 -> R strongly in ^(S^Z^). 
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On the other side, the periodicity of z H> ipj(Sz) in its first variable yields the 
identity 

ttt- / d\Vj dz = (Id+hjG)Sei . 

{S-^p k) \Js-i Z ™ 

By (|7.16[) the left hand side converges to Re±, while the right hand side converges 
to Sei; since \Sei\ =/= 1 = |-Rei|, this is a contradiction. 

The argument for (|7.13[) relies on a bending ansatz for inextensible rods. Let 
va = (vq, Vq, Vq) denote a smooth, (1 + s)-periodic curve from R to R 3 with \v' Q \ = 1 
and 

(7.17) Vo(0) = (0,i,|), «o(l + «) = «o(0) + ei, v' (0) = v' (l + s) = e x . 

Further, we assume that the curve is parallel to the (x\, X2)-plane, i.e. Vq = f . Let 
to := ffj denote the tangent, n :— —t a A e 3 the normal and n :— t ■ n' the signed 
curvature of the curve. Now, define the scaled curve 

v {k) {yi) := kv {yi/k) 

and let t^ k \n^ and denote the associated tangent, normal and curvature. 
Note that the curvature scales as 

(7.18) «<*5(vi) = |«o(l/i/fc). 

Based on we construct a 3d deformation. For y <E Zp set 

4%) : = „(*)((! + S ) yi ) + ( y2 _ ±)n«((l + S ) yi ) + (2/3 - |)e 3 + MZe lVl . 

By construction the map V'(y) := u ^h\y) ~ (Id +hG)y belongs to Xp k \ Thus, by 
definition of m p (Id +hG) we have 

(7.19) mfXU+hG) <\f W (Vu{ k) (y)S) dy. 

A direct computation shows that 

V^ fe) = [t^\n^\e 3 ^S- 1 + (y 2 - §)( {n^)'®e x )S~ X + hG{e x ®ei), 

Multiplication with the rotation R^ :~ (t^ \ \ e^) and a short calculation 
yields the identity 

(R^fVu^S = (Id + i^W + hBW) 

with A^(y) := k^ k) { yi ){y 2 - |) (ei®ei) and B<W := (R^^G (e^e^S. Since 
1 2/2 — 1 1 < p and by appealing to (|7.18p and the smoothness of Vq, we find that 

max \A [k) {y)\<c and \B {k) \ < c\G ei \ 

yez { p k > 

where c only depends on the radius p and the curve Vq. Hence, for small h and 
large k the matrix (i?( fe )) T Vzi^ fc ' ) (y)S is uniformly close to Id and we can appeal to 
condition (W4) (and the frame indifference of Wo) to justify the following expansion: 
For h < h and h sufficiently small there holds (uniformly in y 6 zf^) 

W (Vu {k) S) = Wo((R {k) ) T Vu {k) S) = Woild+k- 1 ^ + fcflW) 
= QoOrV fe) + hB<- k) ) + o{h 2 ) 



28 



S. MULLER & S. NEUKAMM 



where Qo denotes the quadratic form associated with the expansion of Wo at iden- 
tity. Since j B (1/2 — 5) ^(2/2, 2/3) = by symmetry the components of and B^ 
are orthogonal in L 2 (B p ) and we get 

Qoik-'A^ +hB^)d(y 2 ,y 3 ) = ± f Q Q (A™) d(y 2 ,y 3 ) + h 2 \B p \Q a (B^). 

J Bp 

In conclusion, we have 
hmsupi/7 W (Vu^S)dy< ^ if Gei = 0, 



IToo^ fc 7o 7i3 p " uv " ft ~" [2\B p \h 2 max B Qo(B) + o{h 2 ) else, 

where the maximum is taken over the set of all cluster points of the bounded 
sequence (B^) C M 3 . In combination with (f77T9|) this implies (f7TT3|) . 

gtep 3. (Properties of I ( k) ). Set 

m^(F)~ inf ■ A k \ip;F). 

LetQ :={G = -g(e 2 ®e 2 ) ■ g G (0, 1) }. We claim that 

(7.20) VG g M 3 , fc G N : lim sup ^m k) (Id +/1G) < ^Qo(G), 

h-K) h 2 

(7.21) VG G G, k G N : lim inf ^ ( Id > 

/40 h 2 

m ( k ) (Id 4-/lG) 

(7.22) VG E Q : lim sup lim sup V = 0. 

HO fc->oo h 

Argument for (|7.20|) : By definition of m Q k ^ and condition (W4) we have 

im^Id+M?) < J^J^ k) ^o(H+^G) da; = Iq (G) + o(h 2 ) 

and (|7.20l) follows by passing to the lim sup h on both sides. 

Proof of (|7.21[) by contradiction: Suppose that there exist k G N, G G Q and a 
vanishing sequence of positive numbers (hj)j^ such that 

(7.23) lim -Km k) (Id +hjG) = 0. 

.?->■ 00 /ij 

By definition of mj*' there exists a sequence (ipj) C -Xq such that 
lim / W (Id+^(G + V^-)) dy = 0. 

By appealing to the non-degeneracy condition (W3) and geometric rigidity (see 
Theorem 13. ip there exist rotations Rj such that 

(7.24) lim h^Wld+hAG + V^j)- R]\\ T 2, 

j—too \ I 
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Since ipj is periodic in its first and second variable, we get for i = 1, 2 
hfKld-R^eil <K X | (Id +hjG - Rj)*] + \Ge. t \ 



iZ^r 1 [ (M+hjiG + Vri-Rjtedx 
Jzx ' 



IGe, 



<hf\\ (Id +hj(G + Vipj) - Ri)ei 1 1 L2 (z w j|^* } r 1/a + I Ge, | 

Hence, in virtue of (I7.24[) . this implies | (Id. — Rj)&i\ < chj for % = 1,2 and a con- 
stant c independent of j. By appealing to the property that Rj £ SO (3), we find 
that | Id. — -Rj j < 2chj. Consequently, by (|7.24[) the sequence (V<fj) is bounded in 
L 2 (Z^ k) ). Now, we can proceed as in the proof of Theorem ll.il (Step 2) and get 

1 1 f 

liminf -jm^Id+^G) > inf -r / Q (sym(G + Vcp)) dy. 

(k) 

By appealing to the periodicity properties of maps in Xq and since Qq is positive 
definite on the subspace of symmetric matrices, we can easily derive a (non-optimal) 
lower bound; namely, the right hand side is bounded from below by 

c^|G : (e 2 <g>e 2 )| 2 for some positive constant c (fc) > 0. 

For G £ Q this expression is strictly positive — in contradiction to ()7.23|) : thus, 
([731]) follows. 

The argument for (|7.22|) is more subtle and relies on the observation that Iq 
can be written as the clastic energy of a membrane occupying the slender domain 

OW :=wx(0,i), w:=(0,l) 2 . 

To make this precise, define the functional 

and notice that with the scaling ii(x) = ku(x/k) we have 

j^( u ) = ^l^w (m)d X . 

As a consequence it follows that 

(7-25) ^ ) ^ +/ ^ <inf{^)(,): u £ X%> +hG } , 

where X [ p ] := j u £ W 1,2 (fi (fe) ; R 3 ), u{z) = Fz on du)x(0, i) | . 

In [T3] it is shown that T-converges as k — > oo to a membrane energy which 
is zero for short maps, i.e. for deformations u : uj — > R 3 with V« T V« < Id. Since 
(Id +hG) T (Id +hG) < Id for all G £ Q and sufficiently small h > 0, it is natural 
to prove (|7.22|) by appealing to a suitable recovery sequence. For our purpose it 
is not necessary to derive the precise form of the limiting functional. Indeed, the 
following is sufficient: For all G £ Q and 0<ftcl there exists a sequence (u^), 



u (k) g X$ +hG such that 



limsup J (fe) (u (fe) ) = 0. 

k— >oo 
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In combination with (|7.25|) this clearly proves the assertion. The existence of such 
a sequence is shown for instance in |7| 114] . 

Step 4- (Conclusion). Argument for a). We only prove (|W3'I) . Since Wq is non- 
negative and non-degenerate (see condition (W2)), we have for all ip £ Wp£(Y; K 3 ) 

/ W{y,F + VV>) dy> f W (F + W) dy > c' j dist 2 (F + Vtp, 50(3)) dy. 

JY J Yq J Yq 

Here and below, c',c",c"' are positive constants that only depend on Wo and on 
the geometry of Yq. By geometric rigidity (see Theorem 13. ip we get (for some 
R £ 50(3)) 

/ W{y,F+Wtp) dy > c' f \F + Vt/j-R\ 2 dy > c"\F-R\ 2 > c"'dist 2 (F ) 50(3)). 
Jy Jy 

Inequality (★) is valid, since 

inf / \G + V?A(j/)| 2 dy > 
g ^Jy 

where the infimum is taken over all Gel 3 with \G\ = 1 and ip £ W^ 2 (Y;R 3 ) 
satisfying J Q f Q d 3 ip dy 2 dyi = 0. 

Proof of b) by contradiction. Suppose that there exists F £ M 3 and k £ N such 
that W$ m (F) = 0. Then (fTTTT]) implies that m {k) (F) = m {k) {F) = 0. As in the 
proof of ([7T3| . we find that m Q k) (F) = implies \Fe x \ = 1, while m ( p k) (F) = 
implies |-Fei| > 1. But this is a contradiction. 

Statement c) is a direct consequence of (|7.13|) . the splitting in Step 1 and the 
property that W{y, Id) = Wo (Id) = for y £ Y Q . 

Statement d) follows from (|7.13p . the splitting in Step 1 and (17.201) . 

Proof of e). Let G £ Q, i.e. G = —g{e2®e2) where g £ (0, 1). Note that by 
Step 1 we have W^^Id) = m Q k \ld) + (Id) = mj^ (Id) and therefore 

(7-26) <i(Id +hG) - (Id) - m[ k) (Id +hG). 

With (TT2"0)) and ([7T2T|) we obtain 

n r ■ f m k) (ld+hG) _ m {k) (Id +hG) ^ 1_ 

< hmmf u \ 9 < hmsup u v , 9 < -Qo{G). 

HO h z hi Q h z 2 

Hence, if G M 3 fulfills JTJSJ it must satisfy (tr^ , (?) = and we get 

" 1? ( W ^ U +hG) ~ ^(W) - h (a« , G)) > 0. 
On the other side, the combination of (17.26[) and (|7.22[) implies that 

^W^(ld + hG) = 0, 
which completes the proof. □ 
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